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ABSTRACT 

An examination of a mathematical analysis of Pulse Tube refrig- 
eration as proposed by Rea was carried out. Rea originally investi- 
gated the performance of thermal regenerators subjected to rapid 
pressure and flow cycling and found that this class of regenerators 
is characterized by two basic parameters; 

(1) The ratio of the convective heat transfer within the regenerator 
to the heat generated by compression, 

(2) The ratio of the rate of mass storage within the void volume of 
the regenerator to the through-flow. 

Rea's proposal was to treat the Pulse Tube as a very efficient regen- 
erator in series with a very inefficient regenerator, and, thus, 
characterize the performance of the Pulse Tube by the same two para- 
meters. 

The examination of Rea's proposal indicates that, in addition 
to the two parameters which describe the performance of an efficient 
regenerator, there is, for the case of the inefficient regenerator, a 
third parameter which depends upon the average pressure and the rate 
at which the pressure and mass flow cycle. 

An experimental evaluation of Rea's proposal modified to include 
the third parameter was performed on a Pulse Tube consisting of a 
spherical matrix regenerator in series with a void tube. There was 
good agreement observed between theoretical and experimental results* 



Thesis Supervisor? Joseph L. Smith, Jr. 

Title: Associate Professor of Mechanical Engineering 



-ii- 



TABLE OF CONTENTS 



Pag© 

ABSTRACT ii 

LIST OF TABLES v 

LIST OF FIGURES vi . 

ACKNOWLEDGEMENTS viii 

NOMENCLATURE ±x 

CHAPTER 

I. INTRODUCTION 1 

General , 1 

Previous Work 3 

II. EXAMINATION OF REA’S THEORY AS APPLIED TO A REGENERATOR 

AND "pulse Tube" 6 

Rea's 2-Part Model 6 

Basic Differential Equations for the Regenerator 

and "pulse tube" .......... 7 

Set-up of the Differential Equations for 

Computer Solution 13 

Method of Solution 15 

III. ANALYSIS OF THE PULSE TUBE 17 

The System 17 

Analysis 18 

IV. EXPERIMENTAL PROCEDURE 23 

Description of Apparatus .............. 23 

Method of Testing ..... .... 26 

Sample Calculation 3h 

V. RESULTS 36 

Presentation of Results 36 

Discussion of Results U 6 

VI. CONCLUSIONS AND RECOMMENDATIONS k9 



-iii 



TABLE OF CONTENTS (continued) 



Page 

Conclusions 1 j9 

Recommendations 5>0 

REFERENCES 52 

APPENDICES 

A. REA'S 2-PART MODEL AND THE DEVELOPMENT OF THE GOVERNING 

DIFFERENTIAL EQUATIONS Sk 

General 5U 

SECTION 

Aol REA'S 2-PART MODEL AND DEVELOPMENT OF THE REGENERA- 
TOR DIFFERENTIAL EQUATIONS 55 

Rea's 2-Part Model 55 

Development of the Governing Differential 

Equations for the Regenerator 56 

Simplification of the Differential 

Equations 62 

A. 2 THE GOVERNING DIFFERENTIAL EQUATIONS FOR REA'S 

2-PART MODEL AS APPLIED TO A "pulse tube" 7h 

A. 3 METHOD OF SOLUTION FOR THE DIFFERENTIAL 

EQUATIONS 8U 

Regenerator 81* 

"pulse tube" 86 

B. COMPUTER PROGRAM FOR THE SOLUTION OF THE "pulse tube" 

EQUATIONS 88 

General 88 

Computer Program 89 

C. THEORETICAL AXIAL TEMPERATURE AND AXIAL MASS FLOW 

DISTRIBUTIONS, 2 -PART MODEL 90 

D. DESCRIPTION OF THE TEST APPARATUS 10i| 



-iv- 



TABLE 



LIST OF TABLES 



Page 

I. EXPERIMENTAL AND THEORETICAL RESULTS - SERIES 1 . 37 

II. EXPERIMENTAL AND THEORETICAL RESULTS - SERIES 2 38 

III. EXPERIMENTAL AND THEORETICAL RESULTS - SERIES 3 ...... . 39 



-V' 



LIST OF FIGURES 

FIGURE Page 

I. SCHEMATIC OF THE PULSE TUBE 1 

II. SCHEMATIC OF REA'S 2 -PART MODEL FOR A REGENERATOR ..... 1* 

III. WORKING SCHEMATIC OF THE PULSE TUBE 17 

IV. SCHEMATIC OF TEST APPARATUS 21* 

V. TYPICAL "pulse tube" PRESSURE- TIME HISTORIES 32 

VI. SERIES 1 - RESULTS 1*0 

VII. SERIES 2 - RESULTS la 

VIII. SERIES 3 - RESULTS Ii2 

IX. "pulse tube" TEMPERATURE DISTRIBUTION 1*3 

X. "pulse tube" TEMPERATURE DISTRIBUTIONS 1*1* 

XI. CROSS CURVES, C ^ versus <p , m versus § , for 

0 I 0.77 and w z £.0 1*5 

XII. WORKING SCHEMATIC OF THE REGENERATOR 56 

XIII. SCHEMATIC OF THE SYSTEM TO BE ANALYZED FOR THE 

"pulse tube" 75 

XIV. ILLUSTRATION OF COMPUTER SOLUTIONS FOR REGENERATOR .... 85 

XV. ILLUSTRATION OF COMPUTER SOLUTIONS FOR "pulse tube" .... 87 

XVI. THEORETICAL AXIAL TEMPERATURE DISTRIBUTIONS, 2-PART 

MODEL, h = Km 0e ^ 9 , X = 1.67 91 

XVII. THEORETICAL AXIAL MASS FLOW DISTRIBUTIONS, 2-PART MODEL, 

h = Km 0e ^ 9 , X r 1.67 92 

XVIII. THEORETICAL AXIAL TEMPERATURE AND AXIAL MASS FLOW DISTRI- 
BUTIONS, 2-PART MODEL, §10, h-*K, X = 1.67 ... . 93 

XIX. THEORETICAL AXIAL TEMPERATURE AND AXIAL MASS FLOW DISTRI- 
BUTIONS, 2-PART MODEL, § Z 1.0, h = K, * = 1.67 . . . . 9h 

-vi- 



LIST OF FIGURES (continued) 



FIGURE 



Page 



XX. THEORETICAL AXIAL TEMPERATURE AND AXIAL MASS FLOW DISTRI- 
BUTIONS, 2-PART MODEL, (!) = 2.0, h Z K, X = 1.6? . . 

XXI. THEORETICAL AXIAL TEMPERATURE AND AXIAL MASS FLOW DISTRI- 
BUTIONS, 2-PART MODEL, <£ = 3.0, h = K, X = 1.67 . . 

XXII. THEORETICAL AXIAL TEMPERATURE AND AXIAL MASS FLOW DISTRI- 
BUTIONS, 2-PART MODEL, <£ = U.O, h = K, * = 1.67 . . 

XXm. THEORETICAL AXIAL TEMPERATURE AND AXIAL MASS FLOW DISTRI- 
BUTIONS, 2-PART MODEL, c£ = 5.0, h - K, X = 1.67 . . 

XXIV. THEORETICAL AXIAL TEMPERATURE AND AXIAL MSS FLOW DISTRI- 
BUTIONS, 2-PART MODEL, § = 6.0, h = K, X = 1.67 . . 

XXV. THEORETICAL AXIAL TEMPERATURE AND AXIAL MASS FLOW DISTRI- 
BUTIONS, 2-PART MODEL, J = 7.0, h = K, Jj" = 1.67 . . 

XXVI. THEORETICAL AXIAL TEMPERATURE AND AXIAL MASS FLOW DISTRI- 
BUTIONS, 2-PART MODEL, = 8.0, h = K, X = 1.67 . . 



95 

96 

97 

98 

99 
100 



101 



XXVII. THEORETICAL AXIAL TEMPERATURE AND AXIAL MASS FLOW DISTRI- 
BUTIONS, 2-PART MODEL, ^ = 9.0, h = K, X = 1.67 .... 102 

XXVIII. THEORETICAL AXIAL TEMPERATURE AND AXIAL MASS FLOW DISTRI- 
BUTIONS, 2-PART MODEL, $ = 10.0, h s K, X = 1.67 ... . 103 



XXIX. INLET/EXHAUST MANIFOLD 



105 



XXX. SCALE DRAWING OF PULSE TUBE 



106 



XXXI. HEAT EXCHANGER 



109 



-vii- 



ACKNOWLEDGEMENTS 



We wish to thank Professor Joseph L. Smith, Jr. for his 
continual advice during this project. 

We also wish to thank Samuel N. Rea for his many lively 
discussions on the subject at the start of the work. 

This work was supported by and performed at the Cryogenic 
Engineering Laboratory at M.I.T., whose staff members were always 
available for advice and assistance. 

The computer work associated with this project was done at 
the Computation Center of M.I.T. 



NOMENCLATURE 



cross-sectional area of the regenerator or pulse tube, 
ft 2 . 

cross-sectional area of the regenerator or pulse tube open 
to flow (A 0 = £ A), ft 2 . 

heat transfer surface area per foot of length of regenerator 
or pulse tube, ft 2 /ft. 

dimensionless variable defined by equation (A »U8 ) « 
dimensionless constant defined by equation (2. 31). 
integration constant (C 0 a s 

constant defined by equation (2.19), lb.. - °F/hr. 
specific heat of regenerator matrix or pulse tube wall, 

BTU/lbjjj - °F . 

constant pressure specific heat of the gas, BTU/lb m - °F* 
constant volume specific heat of the gas, BTU/lb m - °F. 
gas enthalpy at a position x along the regenerator or pulse 
tube, 3TU/lb m . 

heat transfer coefficient at position x along the regenera- 
tor or pulse tube, BTU/hr-ft 2 - °F . 
heat transfer coefficient, BTU/hr-ft 2 - °F. 
thermal conductivity, BTU/hr-ft-°F« 
regenerator or pulse tube length, ft. 
gas mass flow at position x along regenerator or pulse 
tube, lb^/hr. 

dimonsionless mass flow (m s m/m 0 ). 
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o 




n 



p 

Q 

R 

t 

T 



T 

g 



u 

V 

V H 

V 

V R 

v T 



X 



mass flow at x r 0 in regenerator or pulse tube, lb m /hr. 
matrix mass per foot of regenerator length or mass per 
foot of pulse tube length, lb m /ft. 
positive exponent whose value depends on regenerator 
matrix (n = 0 for laminar flow pulse tube). 

r\ 

gas pressure in regenerator or pulse tube, lbf/ft • 

heat flux, BTU/hr. 

gas constant, f t-lb^./lb m - 0 F . 

time, hours. 

temperature of the regenerator matrix or pulse tube wall 
at position x, °F. 

temperature of tho gas at position x along the regenera- 
tor or pulse tube, °F. 
internal energy of the gas, BTU/lb m . 
volume, ft^. 

dead volume at end of pulse tube, ft^. 

volume of pulse tube, ft^« 

void volume of regenerator, ft^. 

total volume of tho regenerator or pulse tube, ft^. 
dimensionless distance along regenerator or pulse tube 
defined by equation (2.26). 
distance along regenerator or pulse tube, ft. 



Greek Letters 

<X - dimensionless Pulse Tube parameter defined by equation (3.18). 

P - hoat transfer area per unit total volume of the regenerator 

or pulse tube ( ^ ~ ~i), ft^/ft^© 



-x- 



€ 


dimensionless regenerator or pulse tube void fraction. 


1 1 


dimensionless ratio of specific heats ( X = c /c ). 

.H k * 

dimensionless pulse tube parameter defined by equation 


9 


(2.35). 

gas density, Ib^/ft^. 


e 


dimensionless matrix or pulse tube wall temperature 

(© = t/t 0 ). 


r 


duration of cycle, hours. 


i i 

v-° v-° 


duration of a compression period, hours, 
duration of an expansion period, hours. 




percentage of the total cycle time devoted to compression 
or expansion. 



Subscripts 



avg 


average . 


c 


compression. 


C 


conditions at cold end of Pulse Tube 0 


cond 


conduction. 


e 


expansion. 


eff 


effective . 


exper - 


experimental. 


g 


gas. 


heater - 


0 

used with Q to indicate heater input to Pulse Tube. 


H 


conditions at hot end of pulse tube. 


losses - 


• 

used with Q to indicate conduction and radiation losses 
in Pulse Tube. 


m 


refers to regenerator matrix or pulse tube wall. 



-xi- 



max 


- 


refers to maximum value of the variable 


w or p. 


min 


- 


refers to minimum value of the variable 


P. 


0 


- 


pulse tube or regenerator conditions at 


o 

11 

£ 

II 

X 


pt 


- 


pulse tube. 




rad 


- 


radiation. 




R 


mm 


o 

regenerator; also refers to refrigeration (Q^). 


the or 


- 


theoretical. 




W 




conditions at warm end of regenerator. 
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CHAPTER I 



INTRODU: 



G:-nsrd 

A recent paper by Gifford and Longsworth (l)’" describes an inter- 
esting thermodynamic effect which frequently is ignored in pracxice, 
tut which has a potentially useful application,. As Gifford and 



FIGURE ± 



SCHEMATIC OF THE PULSE TUBE 





Longsworth indicate (2), it is not generally recognized that the cycling 
of a gas in a closed volume can establish large temperature gradients 
within such a volume. Gifford and Longsworth, by suitably arranging a 
u Normal regenerator, heat exchangers, and a void volume, have created a 



v Numbers in ( ) refer to References. 



device which preserves these temperature gradients in an essentially 
static state, even though the gas is undergoing rapid cycling* They 
have celled this device a Pulse Tube. It is illustrated schematically 



in Figure I. 

The construction and operation of the device are simple. This 
simplicity when coupled with the promising experimental results obtained 
thus far, indicates that the Pulse Tube should find practical applica- 
tion in the near future, perhaps as a small cooler in electronic equip- 
ment. C-ifford and Longsworth (2) report refrigeration across temperature 
intervals as large as 2G9° F. with a single-stage system. The units may 
also be operated in multi-stages, so larger temperature differences are 
entirely feasible. 

The Pulse Tube consists of a regenerator connected in series with 
a void tube, which is the pulse tube.'"''’ The void volume of the regen- 
erator is The volume of the pulse tube is V , . Gas enters and 

exits the device through the warm end of the regenerator, which is at 
temperature The volume which is connected to the end of the 

pulse tube, is maintained at a constant temperature, Tp. This may be 

accomplished by circulating cooling water around the volume V_ T . 

n 

The device operates in the following manner: the inlet valve 

opens admitting high pressure gas into the regenerator, thereby 
inducing a mass flow into the regenerator and pulse tube. 



Pulse Tube capitalized refers to the entire device (regenerator 
plus pulse tube); pulse tube without capital letters refers to 
the void tube. Wherever the void tube is intended, but gram- 
matically capital letters are required, the following notation 
will be used: "pulse tube". 
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Concurrently, there is a continuous pressure rise throughout the Pulse 
lube o At so:ze predetermined maximum nressure, p , within the device, 
the inlet valve closes and the exit valve opens. This allows the gas 
within the Pulse Tube to flow back out through the warm end of the regen- 
erator. 

When the inlet valve is open, the gas in the pulse tube is heated 
as it is compressed by the entering gas. Some of this heated gas flows 
into the volume V , where it is cooled by the cooling medium to the tem- 
perature T^. When the exit valve opens, the gas is cooled further by 
expansion as it flows back out of the volume V . Some of this cool gas 
must precool the entering warm gas during the next compression half- 
cycle. The remainder of the cool gas is available for refrigeration. 
Refrigeration becomes available, once a steady cyclic state is reached, 
at the connecting tubing between the regenerator and the pulse tube 0 
This refrigeration load is labeled Q in Figure I and is available at 
the temperature T^. 

Previous Work 

The first analysis of the pulse tube was given by Gifford and 
Longsworth in their first and second papers (1,2). This analysis did 
not predict compatible pulse tube capacity results with those they had 
obtained experimentally. Gifford and Longsworth recognized that their 
original analysis of the pulse tube was inadequate. Their original 
analysis neglected any heat exchange between the gas and the pulse tube 
wall. They discovered that this heat exchange could contribute a sub- 
stantial amount to the refrigeration effect. In other words the pulse 
tube itself acts as an inefficient regenerator. The third paper by 
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Gifford and Longsworth (3) discussed the aspects of this heat exchange 
between the gas and the pulse tube wall in detail. However, no complete 
analysis including this effect was presented. The problem was primarily 
one of determining the appropriate temperature difference to use in 

FIGURE II 

SCHEMATIC 0? REA’S 2-PART MODEL FOR A REGENERATOR (U) 




order to calculate correctly the heat transferred in the process under 
consideration. 

Concurrently with the work of Gifford and Longsworth, Rea was 
studying the performance of thermal regenerators subjected to rapid 
pressure and flow cycling (U). Rea’s objective was to develop a theory 
for predicting the axial temperature distribution, mass flow rate 



distribution, and effectiveness of thermal regenerators operated under 
the above conditions. In developing his theory Rea established his so- 
called 2-Part Model (ij), which is patterned after the Gifford-McMahon 
cycle (5>). A schematic of Rea's 2-Part Model is shown in Figure II* 

By making energy and mass balances for the gas and the matrix in 
an elemental control volume, Rea was able to •write the governing differ- 
ential equations for his problem. After applying appropriate simplifi- 
cations, Rea reduced his problem to a system of two differential equa- 
tions, which he solved numerically with the aid of a computer. A complete 
description of Rea’s 2-Part Model and the development of the governing 
differential equations are presented in Appendix A. 

"When the work of Gifford and Longsworth became known to him, Rea 
recognized that the Pulse Tube operates in the same manner as that 
described by his 2- Part Model, Rea concluded that since a pulse tube 
was nothing more than an inefficient regenerator, his theory should 
apply equally well to the entire Pulse Tube as to a simple regenerator,. 
Thus, Rea extended his theory to include the pulse tube. 

Rea tested his theory for the regenerator and got fairly accurate 
results. Rea, however, did not test his theory as extended to include 
the pulse tube. Thus, it is the purpose of this thesis to determine 
whether or not Rea’s Theory correctly describes the performance of the 
Pulse Tube by, first, examining his theory and analysis for the entire 
Pulse Tube and by, second, measuring the refrigeration effect of an 
operating Pulse Tube. 
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CHAPTER II 



EXAMINATION OF REA'S THEORY AS APPLIED 
TO A REGENERATOR AND " pulse tube" 

Rea 1 s 2 -Part Model (U) 

When faced with the task of formulating a mathematical analysis 
of a given system, the most important step is the selection of an 
appropriate model « Ideally one wants a model that is amenable to a 
theoretical treatment. On the other hand one does not want his model 
to be over-simplified, for then the problem may lose its significance* 
Rea selected as his model for the regenerator, and, hence, the pulse 
tube, a model that satisfies both of these criteria. This "2-Part 
Model", as Rea refers to it, is shown schematically in Figure II of 
Chapter I. The usefulness of this model is proved by the fairly 
accurate results which are obtained by its use and by the fact that 
certain dimensionless parameters which are definite aids in defining 
the problem evolve from its use. 

The operation of the 2-Part Model is as follows: the inlet 

valve opens allowing high pressure gas to enter at the top of the 
device. During this half-cycle the pressure within the device con- 
tinually increases with the mass flow rate out of the bottom of the 
device equal to m Q , After a predetermined maximum pressure, p 
is reached within the device, the inlet valve closes and the exit 
valve opens. During the exhaust half-cycle the process described 
above is reversed, with mass flow rate into the bottom of the 
device equal to m Q . The pressure within the device decreases to 



■5m 



its minimum value, P ni ^ n , during this half-cycle. After a period of 
time, a steady, cyclic condition is reached. 

In examining Rea’s theory for the Pulse Tube, the theory for the 
individual components, i.e<>, the regenerator and pulse tube, will be 
examined separately. 

Basic Differential Equations for the Regenerator and " pulse tube 1 1 

The governing differential equations are derived in detail in 
Appendix A by making a mass balance and energy balance for the gas and 
for the matrix in an elemental control volume, dx, at a position x 
within the regenerator and pulse tube respectively# The derivation of 
the regenerator equations is entirely attributable to Rea (U). The 
modification to Rea's Theory as applied to a pulse tube is original# 
The basic equations for the regenerator and pulse tube are: 



Gas Energy Balance : 



C P & ♦ “A 



t) + s^o ap 
r at 



1 o , 



( 2 . 1 ) 



Matrix Energy Balance : 



hA(T - T) 
T g 



Mc m d -l 

at 



( 2 . 2 ) 



Gas Continuity : 



a m _ d I P 



— : : - o or - * ° g 



A 0 1 1 p 



ax r at \ t 



R T At t/ a t 
g' ' g g 



( 2 . 3 ) 



■ 7 < 



The assumptions necessarily made in deriving equations (2.1), 



(2.2), and (2.3) are: 

(1) Zero thermal conductivity in the direction of flow for 
both the gas and the matrix and infinite thermal con- 
ductivity in the matrix transverse to the flow direction, 

(2) The pressure drop in the device is neglected, 

(3) The gas obeys the perfect gas law, 

(h) The gas transport properties and the matrix properties 
are constant, 

(5) The cross sectional area of the device open to flow, 
is constant. 

In addition to these assumptions it is necessary to assume that 
the heat transfer area per unit volume of the pulse tube is fairly 
large, in order that the pulse tube can be considered as a regenera- 



Boundary conditions for equations (2.1), (2.2), and (2.3) for 
the regenerator are: 



Boundary conditions for equations (2.1), (2.2), and (2.3) for 
the pulse tube are: 



tor. 




(2.U) 



at x = L: T g = T gW (t) . 




( 2 . 5 ) 



at x = L: T g = TgC (t) 
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Note that these boundary conditions are identical and have been 
presented separately only to show that the position x : L for the 
regenerator is the warm end of the regenerator, while the position 
x = L for the pulse tube is the cold end of the pulse tube, when the 
Pulse Tube is operated as a refrigerator. 

For a regenerator (and pulse tube in this case) operating in the 
steady, cyclic manner described here, certain cyclic conditions must be 
met. At a given position x, these conditions are as follows: 



<£hA T (T g - T)dt Z N Crrl 21 dt = 0 

6 ' t 

<^mdt = 0 



( 2 . 6 ) 

(2.7) 



dt = o 

<^(mT g )dt = C c 



( 2 . 8 ) 

(2.9) 



Equation (2.6) says that over a complete cycle there can be no 
net heat stored by the matrix. Equation (2.7) says that there is no 
net mass stored over a complete cycle. Equation (2.8) states that the 
pressure is cyclic. Equation (2.9) results from equations (2.6) and 
(2.8) applied to a cyclic integral of equation (2.1). C Q is constant 
with both time and position. The cyclic conditions, equations (2.6) 
to (2.9), are general in nature and must be satisfied by any steady, 
cyclic solution to equations (2.1), (2.2), and (2.3). 

The complete cycle will now be broken into two half-cycles of 
equal time duration. Also, the pressure is 'assumed to be a sawtooth 



function of time such that 



dp 

dt 



S a constant. Other approximations and 
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assumptions are as follows: 



(1) The temperature difference between the gas and the matrix 
at any point within the regenerator (or pulse tube) is 

very small. This enables one to substitute T for T in the 

§ 

final continuity relations. 

(2) The product of the average equals the average of the product. 
This approximation is useful in applying the cyclic condi- 
tions, equations (2.6) through (2.9). 

(3) For the regenerator only, 

A /p_WI A 

bt \T g ] T dt 

With these conditions imposed, the final set of equations to be 
solved for the regenerator are: 



•c p (mT ) + hA T (T gc - T) + dp 

dx R dt 



I 0 



( 2 . 10 ) 



c„ 1. (mT ge ) + hA T (T ge - T) - ^o 



dx 



R 



dp 

dt 



= 0 



( 2 . 11 ) 



dm _ A 0 






9 


(2.12) 


dx RT 


dt 






”(Tgc - T 


ge^ 


~ c : 


> 


(2.13) 


T gc * T ge 


= 2T 


0 


(2.11;) 


boundary conditions 


for these equations are: 




at x = 0: 




T - 


T o ’ 








m z 


*0 


(2.15) 


at x : L: 




T = 
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The final set of equations to be solved for the pulse tube are: 



-c. 



J. (mT ) + hA (T - T) + 
dx g 1 gc 



Vo 

R 

C P k ( ” T ge ) * “ T (T ge - T) - 

R 



ge' 



dp 

dt 

dp 

dt 



dm _ ^ 

dx RT 



dp 

dt 



1 - 2C /JLj _E I [Eq 

1 2T— 1 / At T \ m / 



m(T - T ) = C, 
gc ge 1 



+ T__ ; 2T 
gc ge 



r o 



: o 



V 

KA t L/ 



( 2 . 16 ) 

(2.17) 

( 2 . 18 ) 

(2.19) 

( 2 . 20 ) 



The boundary conditions for these equations are: 



at x = 0: T = T 



at x : L: 



m - m„ 



T = T r 



( 2 . 21 ) 



Note that the set of equations for the regenerator is identical 
with the set of equations for the pulse tube with the exception that 
the continuity equations (2.12) and (2.18) are different. The differ- 
ence between these two equations is that, in the case of the regenera- 
tor, the second term of equation (2.3) can be shown to be negligibly 
small in comparison to the first term, whereas in the case of the 
pulse tube, the second term is not small for all conditions of 
operation. Rea assumed that the continuity equation as written for 
the regenerator applied equally to the pulse tube. The details of 
the relative orders of magnitude of the terms of equation (2.3) for 
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the regenerator and pulse tube are presented in Appendix A. 

Equations (2.10) through (2 •lit) are five equations in five un- 
knowns for the regenerator and can be reduced to two equations in two 
unknowns by subtracting (2.11) from (2.10), and then by substituting 
into this result equations (2.12), (2.13), and (2,lli), plus the com- 
monly known relationship R = c p - c v . The result of this manipulation 
is: 

• dT _ hA^ + A 0 

' • — — 

dx 2 cm c„ 

P P 



dp 

dt 



- 0 



( 2 . 22 ) 



and. 



dm A„ dp 

dx ~ RT dt 



(2.23) 



The three boundary conditions (2.15) are required since, although 
in equation (2.22) is a constant, it is not known at this point. 

Idkexiri.se, equations (2.16) through (2.20) are five equations in 
five unknowns for the pulse tube and in a similar fashion to that 
described above can be reduced to two equations in two unknowns. 

These results are: 



• dT hAmC^ A 

m - j 1 + _£ 

dx 2cm c 

P P 



dp 

dt 




2C 






'At T 



f-l/, 





X - 1 



1 



= 0 



(2.2U) 



• 12 - 



(2.25) 



and. 



dm _ A 0 
dx RT 



dp 

dt 



1 - 2C 



l\ P 1 AoV v 

^y-lj At T \m l\ KA t L 



Again, the three boundary conditions (2*21) are required, since C-^ is 
unknown at this point. 

Set-up of the Differential Equations for Computer Solution 

The following dimensionless variables are now introduced: 

Q Z t/t 



m 



- m/mo 



w - 



V 



W 



dp 

dt 



(2.26) 



x 

L 



Also it is assumed that the heat transfer coefficient can be 
correlated as: 



h - Knf 1 

for the regenerator, and as: 
h = K 



(2.27) 



( 2 . 28 ) 



for the pulse tube, which is assumed to be operating in laminar flow. 

When the new variables ( 2 . 26 ), along with the assumption (2.27), 
are introduced into the regenerator equations (2.22) and (2.23), the 
final regenerator equations result in: 



d0 

dw 






m 



2-n 



(2.29) 
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and 



dm _ 1 
dw ~ © 



(2.30) 



where 



c = y-i /ka t l c x 1 



V / 2m„ 



dP 

dt 



(2.31) 



The boundary conditions for the regenerator become: 

at w : 0: m = l 

©= 1 



at w = w, 



max 



V 


dP 




dt 



w 



: 0 : 1:0 



W 



( 2 . 32 ) 



In a similar fashion, when the new variables (2.26), along with 
the assumption (2.28), are introduced into the pulse tube equations 
(2 ,2ii ) and (2.25), the final pulse tube equations result in: 



I 1 * / y-i l 
d 0 _ c 1 e 1 ~ R 1 y / 

dw 



(2.33) 



nr 



and 



dm - 1 



I 



1 -$C 
0m 



dw © 

where C is defined by equation ( 2 ,31 ) and 



(2.3W 






y 

y-i 






KA.^L j At T c 
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(2.35) 



The boundary conditions for the pulse tube become: 



at w s 0: 



at w = w. 



max 



m = 1 

©= 1 

V dp 

T m R dt 
o o 




(2.36) 



Method of Solution 

The method of solution for the regenerator equations (2.29) and 
(2.30) is to treat C as a parameter. This makes the problem an initial 
value problem, since both © and m are known at w = 0. To obtain a set 
of solutions, one must know the ratio of the specific heats, X , and 
the tyoe of matrix used in order that n can be assigned a value. The 
parameter C is then assigned a range of values, and the equations are 
numerically integrated from w = 0 to w = w^. Solutions for n = .59 
and X - I .67 are presented in graphical form in Appendix C« 

The method of solution for the pulse tube equations (2.33) a ^d 
(2.3k) is similar to that for the regenerator. Row, however, there is 
a second parameter, , involved. In this case is assigned a value, 

O - C£>^, and then C is assigned a range of values for this (p . Again 
the eouations are numerically integrated from w z 0 to w = % x . The 
result is a family of curves of © versus w and m versus w, with each 
curve having a specific value of C and each graph representing one 
value of the parameter <£) . Solutions for various values of and 
X z 1.67 are presented in graphical form in Appendix C. 

To use the curves in Appendix C one must know w^^ and Q ^ 
for the case of the regenerator. In the case of the pulse tube, 
w max> ® C* and must be known. Then entering the appropriate 
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graph with w = w raax and © = © VJ , or © = © c , a value of C can be 
read. The axial temperature distribution is determined by following 
that value of C down to w = 0. Also, once C is known, the mass flow 
into the top of the regenerator (or the bottom of the pulse tube, as 
the case may be) is determined by entering the graph of m versus w at 
w - W psy and reading off the value of m corresponding to the inter- 
section of the C curve and the w = w^^ line 0 As in the case of the 
temperature distribution, the axial mass flow distribution is deter- 
mined by following the same C curve from w = w Kax to w = 0. This 
yields the mass flow distribution from the warm end to the cold end 
of the regenerator or from the cold end to the warm end of the pulse 
tube . 

From this point on equations (2.29), (2.30), (2.33), and (2 ,3h) , 
plus the boundary conditions equations (2.32) and (2.36), will be 
referred to as Rea's Modified Theory for the Pulse Tube. 
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CHAPTER III 



ANALYSIS OF THE PULSE TUBE 



The System 

The Pulse Tube consists of a regenerator in series with a pulse 
tube, as shovm schematically in Figure III. This, then, is the system 
which will be analyzed. 



FIGURE IH 

WORKING SCHEMATIC OF THE PULSE TUBE 
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Analysis 



If one considers the control volume around the entire Pulse Tube 
as shown in Figure III and writes an energy balance for one cycle of 
length T', the result is: 



r c 



m T 



m T 
e ge 



W 






Tz o 



(3.1) 



where the subscript W refers to the warm end of the regenerator. In 
equation (3*1) the first terra represents the enthalpy carried by the 
gas into the control volume at the warm end of the regenerator, the 
second term is a known heat quantity introduced into the control volume, 
the third term is the heat associated with conduction and radiation that 
is absorbed by the control volume, the fourth term represents the enthalpy 
carried by the gas out of the control volume at the warm end of the regen- 
erator, and the fifth term is the heat removed from the volume V^. 

Now, if the control volume around the pulse tube only is considered, 
one can write an energy balance for one cycle of length Y as follows: 



*Y c m T I + Q' 'Y - 'Y c 

1 c p \ c gc| r losses e p 



m T 
e ge 



- Q H r r 



(3.2) 



where the subscript C refers to the cold end of the pulse tube and 

• | 

Q , are the radiation and conduction losses associated with the 

losses 

pulse tube 0 

If one now makes the assumption that the period of compression, 
Y Qi is equal to the period of expansion, 'T , and, furthermore, that 
either of these periods is equal to where is the percentage of 
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the total time that is devoted to compression or expansion, then equa- 
tions (3.1) and (3*2) can be ircitten as: 



^ C T 



m T„ - T ' 
gc ge 1 



W 



+ Q' 



+ Q 



heater ^losses 



- Q H = 0 , (3.1a) 



and 



y C r 



m T - T 
gc ge 



+ ^'losses " Q H = 0 * 



(3.2a) 



Ely combining equations (3.1a) and (3.2a) one obtains the result: 



Qr = V C T 



m r 



T - T | 
gc gel c 



" ^ l T gc “ T g el w 



(3.3) 



where 



Q - Q 
R “ heater 






losses 



(3.U) 



and mCp | T^ c - T^ e j is the difference in enthalpy flux at some 
point in the regenerator or pulse tube. This enthalpy flux can be pre- 
dicted by Rea’s Theory as presented in Chapter II . 

The difference in enthalpy flux at the cold end of the pulse tube 
and at the warm end of the regenerator will now be found. 

From equation (2.31), 



C = 2C 
1 





d£ 

dt 



(3.5) 



and from equation ( 2 . 13 ), 



C-[_ = m 





(3.6) 
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Now consider the regenerator* At the bottom, or cold end, x “ 0 and 



w w 

in = m 0 = mQ 0 At the top, or the warm end, x = L, T = ^.j, and m •- m^. 
From equation (3.6), 






T - T 
gc ge 



W 



W 



(3.7) 



Then from equation (3.5)* 



= 2C t 



R 



J-l/ \ka t l 






R 



dp 

dt 



(3.8) 



Now consider the pulse tube. At the connection between the cooler 

• • • 

section and the pulse tube, x = 0 and m = m Q = mjj. At the bottom, 

• • • 

or the cold end, x = L, m = m^ = n^, and T = T^ = T^. From equation (3.6), 

(3.9) 



^lV- T ge] 0 = 1 °X | c = l*, 



Then from equation (3.5)* 



t°x) 



pt 



= 2C 



pt 



x 

y-1. 



\KA t L 



“h 



pt 



dp 

dt 



(3.10) 



By substituting equations (3.7) and (3.9) into equation (3*3) one 
obtains: 



q r z t c P 



( c l)p t - ( c l) 



(3.11) 
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Now by substituting equations (3.8) and (3.10) into equation (3 oil) the 
result is: 



Q r =2^c_ 




d£ 

dt 



V 



^(kAjl) t Cpt " ^ (ka^l j ^ C R 



• ( 3 . 12 ) 



Equation (3«12) is identical with equation (3«U) of Rea (It) with the 
exception of the 4 ^ factor, which apparently was omitted by Rea. The 
method of determining this factor will be commented on in the following 
chapter. 

Equation (3.12) will now be simplified. One is able to write for 
a perfect gas: 



pV = mRT 

o 



(3.13) 



where m is the mass. If this equation is differentiated with respect 
to time for a particular volume, the result is the general relation- 
ship: 



m = V 
R 




P dlj 

V dt 



(3.1a) 



Now consider the particular volume to be V Rj the volume of the cooler 
at the end of the pulse tube. The assumption is now made that there 
is perfect heat transfer in V . This means that the gas in the volume 
V R is at a constant temperature T . With this assumption equation 
( 3 .lU ) becomes: 
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(3.1p) 



* v w 

m H ” 



RT, 



H 



dp 

dt 



From equation (3.l£) it follows that 



* m p V w 
m = C 

c — 



R1 k 



dp 

dt 



(3.16) 



By substituting equations (3*l£) and (3.16) into equation (3.12) and by 
rearranging, equation (3.12) can be written: 



q e = 2cW_L1 2 ± l 

U-i) dt t h Up )pt v 



c pt " ^ C R 



(3.17) 



where 



oc = 





Ipt 


( K , 


Pi 


Ir 



(3.13) 



and 



0 5 



a t lc 



V 



(3.19) 



Equation (3«17) is the governing equation for the Pulse Tube 0 
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CHAPTER IV 



EXPERIMENTAL PROCEDURE 

Description of Apparatus 

Rea* s Modified Theory for the Pulse Tube was tested using an 
apparatus constructed as shown schematically in Figure IV. The regen- 
erator consisted of a 12"-long, stainless-steel tube, 1" O.D. by 
0.020"-wall-thickness, packed with O.O^C'-diameter lead shot. The 
pulse tube consisted of a ^"-long, stainless-steel tube, £/8" O.D. by 
0.0l6"-waH-thickness. The bottom of the regenerator was connected to 
the bottom of the pulse tube by a l/h" O.D. by 0.020"-wall-thickness 
stainless-steel tube. Fine-mesh-bronze screens were installed at the 
bottom of the regenerator and pulse tube to smooth the flow. The heat 
exchanger at the top of the pulse tube was constructed from thin-wall 
1/8 "-copper-nickel tubing. The rate of flow of cooling water around 
these tubes was 0.£ gal/min 0 A heating coil of high resistance wire 
was wound around the connecting tubing between the regenerator and 
pulse tube as a means of introducing a known heat load into the device. 
Copper-constantan thermocouples were attached along the length of the 
regenerator and pulse tube. Pressure transducers were installed at 
the warm end of the regenerator and at the water-cooler end of the 
pulse tube. The entire Pulse Tube was suspended in a vacuum- tight 
chamber for thermal insulation. During data runs this vacuum was 
maintained at approximately 2 x 10”£ mm. Hg. Helium was used as the 
working gas for the Pulse Tube. A four-speed motor allowed a choice 
of cycling rates corresponding to 22. U, Uj*7, 82 . 3 , and l£8.0 
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FIGURE IV 



SCHEMATIC OF TEST APPARATUS 



SUPPLY 



EXHAUST 




VARIAC 



W/-J 
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VACUUM JACKET 



reversals per minute. A complete description of the apparatus is 
contained in Appendix D. 

The important characteristics of the device were: 

(1) V R r 3.27 in3 

(2) V pt = 3.39 in3 

(3) - 0.673 in3 

M v ptA H = 5.0U 

(5) v pApt = 0.96U 

(6) e R r 0.388 

(7) e pt = 1.0 

(8) r 1.67 

(9) (K^) 1 595 k BTU/in-hr-ft 2 -°F 

(10) (K^) r = 1.1U x 105 (n^)*^ BTU/in-hr-ft 2 -°F 

(K^)p*, assuming fully developed laminar flow, was obtained from the 
expression: 

(k/3) - IT^iK 

d h 

where Nu = U .36 and j3 - . 

V 



(K£) was obtained as follows: 



for a spherical matrix. 



R - 6 ( 1-6 ) 
P R " d * 



where d is the diameter of the matrix spheres. 
Rea's experimental results (U). 



K was determined from 
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Method of Testing 



Basically the objective of the experimental work was to obtain 
data from a Pulse Tube that could be compared with the refrigeration 
effect which Rea's Modified Theory predicts for that Pulse Tube. 

One of the most important experimental results was the axial 
temperature distribution, particularly of the pulse tube. From the 
equation (2.26), w is defined as: 



w = 



T m R 
0 0 



dt 



x 

L 



For the pulse tube, from equation (3*15): 



• • Vtt 

v , H = JL 



dp 

dt 



Oul) 



0w2) 



Therefore, for the pulse tube: 



- V x 

’ v H 




(U.3) 



where x is the distance along the pulse tube measured from the top as 
shown in Figure III. For a given pulse tube, equation (U.3) allows 
one to specify locations along its length as a function of Wp^. This 
was helpful in that thermocouples could be located at any distance x 
along the pulse tube and that position could be translated easily into 
a corresponding w^. 

The assumption was made that the temperature of the gas inside the 
cooler was uniform and at the same value as that of the cooling water. 

In reality this is not possible, but by using thin-wall cooler tubes 
and by allowing a short dwell time between the compression and expansion 
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parts of the cycle, "perfect" heat transfer could be approached. In 
this case (^o)pt became simply Then by measuring the temperature 

of the pulse tube at each thermocouple location and by dividing each of 
these temperatures by to convert to the non-dimensional temperature 
variable 0 , the experimental temperature distribution of the pulse 
tube could be determined readily. When this experimental distribution 
was compared with theoretical temperature distributions, one obtained 
typically: 




The shape of all the experimental curves was similar to the shape of 
the theoretical curves except in the vicinity of w^^, indicating that 
the laminar flow assumption for the pulse tube was fairly accurate. 

In normal application the theoretical curves would be entered knowing 
a value of w^^ and ©q. However, in this case, since the entire pulse 
tube temperature distribution was known, it was more appropriate to use 
(e c } eff as an entering argument. This gave a more accurate value of 
C ^ than if © c had been used. 

The selection of the method of comparing the calculated results 



- 27 - 



with the experimental results will now be discussed. From equation 

e 

(3 eh) it should be noted that Q R includes the input power from the 
heater plus the conduction and radiation losses in the regenerator and 
pulse tube* Since both conduction and radiation losses depend upon the 
axial temperature distribution, these losses would be approximately con- 
stant for any operating conditions of the device which resulted in the 
same temperature distribution* This suggested a possible procedure for 
testing wherein a change in the measured value of Q R (Qh ea ^ er flosses) 
could be compared with a change in the calculated value of Q_ for simi- 

JrC 

lar temperature distributions and for identical changes in d£ 

dt 

The Q^osses* course > hod to he determined. As a first approxima- 
tion they were analytically obtained from equation (3*17) for no 
heater input* This approximation assumed that the analytical express- 
ion is correct o This assumption may not be exact, but should suffice 
as a first approximation. The comparisons suggested above could be 



made on a plot of Q R versus 



dP 

dt 



as shown below: 



Or 
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A smooth curve cannot be legitimately drawn through the measured points 
of the versus plot, since there is a third variable, , 

which was not fixed. Another method of conducting the experiment would 

have been to keep constant, and then compare changes in experimental 

• • 

Qr with changes in predicted Q^. 

Eased on the above discussion an experimental procedure was 
selected as follows: First, a cycling speed was arbitrarily selected 

from the four speeds available. Second, a ©^ was arbitrarily chosen. 
The selection of © had one limitation, that being the value selected 
had to be greater than the lowest possible obtainable with the 
device under the limiting conditions of operation selected. This limi- 
tation was necessary so that a range of operation was available wherein 
the same ©^ could be maintained with an increasing heat load. Third, 

a ratio of p„_,/p . was selected which would maintain the chosen © . 

max' *min C 

Initially there was no heater input. Under these conditions the 
refrigeration effect consisted of conduction and radiation losses only. 
It took 60-90 minutes for the temperature distribution to become steady. 



Once steady state was reached, the 

selected e c was noted. 

Next, the heater power and 



dp 

dt 



required to maintain the 



were increased in such a manner 
that the © remained constant. It was found that response time to 

c 

slight changes in pressures was considerably less than the response 
time to changes in the heater input. Therefore, for rapid stabiliza- 
tion, the heater power was increased a fixed amount and then the 
was adjusted to achieve the desired©!-,. Once heater power was intro- 



dt 



duced, it was impossible to maintain the identical regenerator tempera- 
ture distribution, although ©^ did remain nearly constant. This, of 
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course, infers that the Q^osses were not exactly constant for a series 

of data runs. The temperature distribution of the pulse tube on the 

other hand could be very closely duplicated. Once the steady condition 

was achieved, the |d]o| and Qh ea t e r were recorded. By adding the 

^losses thc ^heater* was obtained. ^ continuing this process 

a series of experimental points of Q p versus dp c could be obtained. 

K dt 

In the above procedure Qh ea t e r was determined by measuring 
wattage and presented no difficulty. Because of the strong . dependency 



of the Pulse Tube performance on 



d£ 

dt 



, provisions were made to measure 



the pressure at both the warm end of the regenerator and water-cooler 
end of the Pulse Tube. To check the pressure drop through the Pulse 
Tube, the water-cooler end pressure transducer was first calibrated 
under static conditions using previously calibrated pressure gauges 
which were installed at the warm end of the regenerator as a reference. 
The pressure was then allowed to cycle. Photographs of the oscillo- 
scope traces of the pressure tranducer output were taken for different 
pressure settings on the gauges. From these photographs and the 
pressure transducer calibration curve, the maximum and minimum pressures 
at the water-cooler end of the pulse tube were obtained. For most 
pressure settings the measurements taken at both ends of the Pulse Tube 
were essentially identical. For cycling speeds up to and including 
82.3 reversals per minute, the only measurable pressure drop for the 
data taken was approximately 2 psi at a supply pressure of 1$0 psig. 

For these runs it made no difference where in the Pulse Tube the 
pressure measurements were taken. With the highest speed of l£8.0 
reversals, the pressure drop ranged from 3 psi at £0 psig to 2h psi at 
112 psig. For these runs the p^^ and p were measured on the pulse 
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tube side of the device 



Typical pressure transducer traces are shown in Figure V. From 
the traces it can be seen that the total period can be obtained pre- 
cisely. However, the values for T T , and °l£|are not quite so 

c' e dtl 

obvious o In the analysis of the Pulse Tube it was assumed that 
1 c ~ ^ Q - WT . In keeping with this assumption the following values 

of yJT were decided upon by taking an arithmetic mean of 'T and T • 

c e 



Reversals/tlin 


^^(sec) 




22. 1* 


1.0 


0.372 


Uu7 


o.5 


0.372 


82.3 


0.27 


0.372 


158.0 


o.iU 


0.372 



Because of the difference in slopes of the compression and expan- 
sion parts of the cycle and the slight non-linearity of the exhaust 
history, average values were also used to determine |^jE|. With a simi- 
lar pressure history, Rea (lj. ) reported good results by taking: 



2E 

dt 



Z Pmax “ fynin 



avg 



4t, 



UoS) 



avg 



where 



At *g = rr = i (t . r ) 



(U.6) 



and 



p = maximum cycle pressure 
max 



p^in - minimum cycle pressure, 



A similar relationship for 



& 

dt 



was used in this case. 
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To determine a predicted value for equation (3*17) was used: 



q R = 2 Y[_L 

y-i 



dp 

dt 



2 V, 



H 



C 



r H Upl 



pt 



C pt " ocC R 



x;here oC is defined by equation (3.18). For any Pulse Tube £ £ R , 

and V H would be known, y and T„ would also be known depending on the 



H 



dp 

dt 



were 



gas used and the temperature of the cooling water. and 
determined by the methods previously discussed. was known 

once the pulse tube temperature distribution became known. To deter- 
mine a value of a value of the parameter C§) had to be determined. 
From equation (2.35)* 



<5=2 X X_ E-i- 

^ l Jf-l ] 1 U^/ At T 0 



In the above equation y and j^L_j were known; T 0 corresponds to Tjj 

and was also known; At was taken as 4^T*, which was obtained from 

equation (U*6); p was taken as p , which was defined as: 

avg 



p„„ z p . + Pmax “ Pmin 

^avg ^min 

2 



(U.7) 



Thus (£> was obtained. Now since , ( w max^pt anc * ^^C^eff wer ® a ^ 

known, the appropriate plot of Q versus w for the pulse tube was 

• • 

entered and a Cp^ extracted. This in turn led to a value of mc/rr.jj 

• .« • 

by the procedure described in Chapter II. From m c /m H , m^ was found. 
This led to a value for The next step was to determine (Wjj^)^, 

which was obtained from: 
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(w ) 

max , 



% 



Vr 

V 



H 



m H 



tt.8) 



Equation (U.8 ) follows from the definition of w niax , equation (2.2 6), 
assuming that m^ is constant along the connecting tubing. All quantities 
in equation (1**8 ), then, were known. Now knowing 0^ and ( v niax ) > a 
value of was extracted from the regenerator curves of 0 versus w. 
From equation ( 3 * 18 ) oC was calculated, since all necessary quantities 
were known 0 Finally, was calculated from equation (3*17) • 



Sample Calculation 

Suppose one wants to predict the refrigeration for a Pulse Tube 
with the following specifications: 

= cooling water temperature = E>02°R 

% = 1*77°R 
Vpt^H = 0.673 in 3 
V R = 3.27 in3 



Refrigeration is desired at 376°R. Hence, Tg/Tjj Z 0.75 and %/Tq z 1.27. 
Helium is to be used as the working gas. The cycle rate is to be 82.3 
reversals per minute. The pressures to be used are: 



p max = 60 psi 8 



Pmin = 10 P si e 



p ata = 15 psia 



-3U- 



O is calculated from equation (2.35) to be 2.1*6. From Figures 
XVIII through XXVIII in Appendix C with w max s £.0 and T c /T h s 0.7£, a 
cross curve of versus Cp^ can be obtained. Entering this cross 
curve with (^) z 2.1*6, a value of C ^ - 0.11*7 is read. In an analo- 
gous manner a cross curve of m versus (j 3 can be obtained. From this 
• • 

curve mg/mj j = £.8 is read. Then, from equation (1*.8), 

( w max^ R - 



Now from Figure XVI of Appendix C with w = l.li* and %/T c Z 1.27, 

= 0.85 is read. Suppose also that c< has been calculated to be 
0.00811., then oC z 0*0069« For a cycle rate of 82.3, - 0.27 sec 

and *= 0.372. From equation (1*.5): 

^ “ 185 psi/sec • 

In addition Z lil.6 BTU/in-hr-ft^-°F. 

Then, from equation (3.17): 

Qp£ ■» U ©25 wQ-tts 



\ 
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CHAPTER V 



RESULTS 

Presentation of Results 

Experimental data was obtained using the apparatus and the pro- 
cedure described in the preceeding chapter. Three series of data runs 
were made. In the first series different speeds and widely variant 
back pressures were used for the various runs. For all runs in the 
second and third series a single speed and a nearly constant back 
pressure was used. Tables I, II, and III summarize the results for 
these series of data. Figures VI through VIII illustrate the compari- 
sons between the predicted and the actual refrigeration loads. Observed 
and theoretical pulse tube axial temperature distributions for the three 
series of runs are presented in Figures IX and X. Figure XI shows a 
plot of m and versus ^ for w z (w max )p- t and © = (©c) e ff» 
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EXPERIMENTAL AND THEORETICAL RESULTS - SERIES 1 
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Discussion of Results 



From Tables I, H, and III it should be noted that oc C_ is less 

rc 

than 10$ of C ^ in all cases for this data and for practical purposes 
can be neglected when calculating 0 R . The data indicates that as long 
as the regenerator used in the device is an efficient one, its physical 
dimensions relative to the pulse tube dimensions can be smaller than 
those of the regenerator used in this experiment before its effect would 
have to be considered. If the c4C R term of equation (3.17) is neglected, 
that equation becomes: 



% = 2 V Hh 



or, for a given pulse tube, 

Q r = Constant 



op 

dt 



% 



% M 



Po 



(5.D 



pt 



dp 

dt 



pt 



(5.2) 



where Cp,. is a function of . All the Figures of 0 R versus 



do 

dt 



in 



this chapter are plotted using equation (5.2). As a result, all lines 



dg 

dt 



plots are straight lines pass- 



of $ = constant on the Q R versus 
ing through the origin. 

For a particular series of runs a single temperature distribution 

was used in determining a C for any specific run in that series, 

pt 

since all runs in each series had practically identical pulse tube 
axial temperature distributions. The observed and theoretical pulse 
tube temperature distributions are shown in Figures IX and X. The 
theoretical results for each series are plotted as single curves, 
since the deviation between runs in any one series was less than ?$. 



It is also interesting to note that the experimental curves are in good 



agreement with the theoretical curves, indicating that the laminar flew 
assumption for the pulse tube is a good one. 



From Figures VI through VIII one can determine approximately the 

accuracy of Rea’s Modified Theory. However, the accuracy shown in these 

Figures is no better than the accuracy of the variables involved in 

arriving at these results. From equation (5.2) it can be seen that 

depends on dR 
dt 

upon Q . In determining 



and C . The quantity is, in turn, dependent 



dp 

dt 



and , average values of pressure and 
time were used 0 The methods of determining these average values seem 
reasonable. However, other methods which would have resulted in differ- 
ent average values could also be justified. The importance of these 
results, then, lies not in attempting to show their accuracy, but in 
showing that the trend of the results can be explained by the theory. 

The dimensionless parameters which evolve from the theory are most 
valuable in understanding the problem. For example, consider Runs 2 and 
5 from Series 1 of Figure VI. Note that the 



dp 

dt 



2 

is approximately the 



same for these two runs and yet Q-p for Run 2 is more than twice 

* % exper 



Q. 



■R, 



for Run 5« This can be explained, of course, by the fact that 

exper 

for Run 5 the value of (jj was almost twice as large as ^ for Run 2. 

As a result Cp^ for Run 2 was more than one and one-half times as large 

as Cp^. for Run 5* An indication of how C^_ varies with $ in this case 

is shown in Figure XI. Thus, for Run 5 was approximately one- 

6 

half as large as ^R^eor Run 2* The large difference in <|> for 
these two runs is attributable to the fact that the back pressure for 
Run 5 was four times larger than the back pressure for Run 2; also 
Run 5 had a v|>T of 0.27, while '4"T for Run 2 was 0 . 50 . The effect of 

these quantities on Q can be easily seen by examining equation (2.35). 



Another example of interest are Runs 5 and I 4 of Series 1 , Note 
that the total experimental refrigeration for these two runs is about 



equal, yet 



dp 

dt 



2 



is more than twice as large for Run 5 as for Run Ij. 

The explanation, of course, was that O for Run 5 was almost twice as 
large as ^ for Run lu The values of ^ were different, since t p'T' for 
Run $ was 0.27 opposed to 0.^0 for Run U and the back pressure for Run 5 
was 20 % greater for Rim 5 than for Run 1*. This resulted in for Run U 
being more than one and one-half times as large as Cp^ for Run 5. This 
difference in C .j. due to was enough to counteract the opposite effect 



onQ Rtheor duet ° 



dt 



and, hence, both conditions resulted in pro- 



ducing approximately the same refrigeration effect. 



-I18- 



CHAPTER VI 



CONCLUSIONS AND RECOMMENDATIONS 

Conclusions 

From the analytical model of Chapter II, there appeared three 
principal parameters which characterize the operation of a Pulse Tube 0 
These parameters are: 



c =^i j 



fl_ 

2m 0 



1 

dt 



( 2 , 31 ) 



_ V 

max " ToirioR 



dt 



( 2 . 32 ) 



which were described by Rea (h), plus a third parameter , defined as 



As 2 /JL\ [_JL\ E_ 1 . (2.3?) 

* \x-y pW ** t o 



The parameter (p resulted from the fact that the ratio 



l! k— 

Vg ^)/\i »y 



for the pulse tube was much closer to unity than 



for the regenerator* 



the term 




Thus, the simplifying assumption which neglected 
in the continuity equation (2.3) in the case of 



the regenerator could not be made in the case of the pulse tube. 

From the theoretical and the experimental results obtained from 
the Pulse Tube under consideration, it was apparent that Pea's Kocified 



Theory as presented in Chapter II does describe the performance of a 



typical Pu3.se Tube. There was some deviation between the refrigeration 
effect measured and the refrigeration effect predicted. However, this 
deviation falls within acceptable limits for heat transfer data in the 
cyclic case. 

It was apparent from the data obtained that there was good agree- 
ment between the experimental and the theoretical pulse tube axial temp- 
erature distributions. Ibis leads one to the conclusion that laminar 
flow in the pulse tube is a good assumption. 

From the analytical results obtained from the Pulse Tube, it was 
concluded that the effect of the regenerator on the overall refrigera- 
tion was negligible. Consequently, the regenerator could have been 
made smaller in this case before its effect would have become signifi- 
cant. 

Experimental and theoretical results indicated that there is a 
unique temperature distribution for a pulse tube operating with a 
specified 0 q. 

Recommendations 

There are several areas of Pulse Tube refrigeration which should 
be examined. First, the uniqueness of the axial temperature distribu- 
tions of several pulse tubes with different values of w^^ operating 
with a specific ©q should be investigated over a wide range of <]> , 

If this uniqueness can be demonstrated for a significant range of w^^ 
and § , then a suggested method for presenting the computer solutions 
for the pulse tube differential equations is as follows: 

(1) Plot for a fixed w a series of cross curves of C_. versus 

max P& 

for various values of 0 C . 

(2) Introduce theoretical axial temperature distributions for 



some value of <^) which covers the maximum feasible range in w. 

With the temperature distribution curves described above, the ©^ end 

w max an ^ P u -^ se tube can he used to determine a curve of constant C pf 

By following along this constant curve to the value of Wj,^ for which 

the cross curves of Cp^ versus d) apply, a corresponding value of ©^ 

can be obtained. This value of © p represents a single curve on the plot 

u 

of the cross curves of Cp t and . Thus, the relationship between Cp^ 
and becomes known* 

Second, a method for determining an optimum size regenerator to be 
used with a given pulse tube should be devised. 

Third, an optimizing procedure for the design of a pulse tube using 
Rea's Modified Theory should be developed. 



I 
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APPENDIX 



APPENDIX A 



REA'S 2-PART MODEL AND THE DE7FL0PMENT OF 
THE GOVERNING DIFFERENTIAL EQUATIONS 



General 

As a mathematical model for his analysis of thermal regenerators 
subjected to rapid pressure and flow cycling, Rea devised his 2-Part 
Model* From this model Rea developed the differential equations that 
apply to regenerators subjected to the above conditions* Rea proceeded 
to apply certain assumptions which simplified these equations* The 
final result was a system of two differential equations, which Rea 
solved numerically. He then reasoned that since a pulse tube is 
simply an inefficient regenerator, the same differential equations 
should apply to it. However, not all of the assumptions which led to 
the regenerator differential equations hold true in the case of the 
pulse tube. Thus, the final differential equations which govern the 
pulse tube are not the same as the differential equations which govern 
the regenerator. This Appendix is broken into three sections. In 
Section A.l Rea's 2-Part Model and his development of the differential 
equations for the regenerator will be presented in detail as in (U). 
Section A *2 deals with the development of the differential equations 
for the pulse tube. The latter section will not be presented in as 
rigorous a manner as Rea's development in Section A.l, since the two 
developments are similar. All assumptions made in each case will be 
described clearly. Section A. 3 discusses the method of solution for 
the regenerator and pulse tube differential equations* 



REA’S 2-PART MODEL AND DEVELOPMENT 



OF THE REGENERATOR DIFFERENTIAL EQUATIONS 
Rea’s 2-Part Model 

Rea’s 2-Part Model, which is shown in Figure II of Chapter I, is 
defined as a regenerator in which the mass flow rates at each end of 
the regenerator are in phase. By "in phase" it is meant that during 
each half-cycle the flow is always unidirectional within the regenera- 
tor, i.e. during the pressurization half-cycle, mass flows in at the 
top of regenerator and exits at the bottom, and during the expansion 
half-cycle, mass flows in at the bottom and exits at the top. 

The operation of Rea’s 2-Part Model is as follows: the inlet 

valve opens admitting high pressure gas into the regenerator. The 

pressure increases within the regenerator and the mass flow rate, m o , 

«* 

exits at the bottom. The mass flow rate at the top, m^j, although 

o 

unknown, is certainly greater than m Q since some gas must enter and 

occupy the volume within the regenerator, then some maximum pressure, 

? , is reached, the inlet valve closes and the exit valve opens, 

max* * 

During the expansion half-cycle mass flow equal to m 0 enters at the 

♦ 

bottom and mass flow equal to m^. exits at the top. During this half- 
cycle the pressure decreases to its minimum value, p^_.^. Once a 
steady, cyclic state is reached, the above processes are repeated con- 
tinuously x-/ith the temperature at the top equal to T^. and the tempera- 
ture at the bottom equal to T q . 

It is interesting to note that, since ir^ is always greater than 

u , all the gas entering at the top does not reach the bottom. Thus, 
Cr 



it is possible for some gas to remain always within the regenerator 
and flow back and forth as the pressure cycles. 



FIGURE XII 

WORKING SCHEMATIC OF THE REGENERATOR 




Development of the Governing Differential Ebuations for the Regenerator 

The development which follows with regard to the regenerator 
differential equations is entirely attributable to Rea (lj). 

Rea's initial assumptions made concerning the development of the 
regenerator differential equations are as follows: 



(1) 

( 2 ) 
( 3 ) 



The pressure drop in the regenerator is neglected. 



The gas obeys the perfect gas law. 

The gas transport properties and the matrix properties are 
constant, 

(ii) The cross sectional area of the regenerator open to flow, A , 
is constant* 

Figure XII shows a schematic of the system to be analyzed* The 

system consists of a regenerator of length L, which is filled with a 

matrix material of large heat capacity, e.g* spherical lead shot. The 

subscript m is used to denote the matrix,. The void volume within the 

regenerator is denoted by V and is equal to A L, Curing each cycle 

of operation the gas in the regenerator is alternately compressed and 

expanded between the pressures p and p . . 

max min 

In accordance with the assumptions listed previously, energy equa- 
tions will be written for the gas and for the matrix contained within 
the control volume of length dx, as shown in Figure XII« 

Gas Energy Equation ; 

Energy flux into _ Energy flux out of Rate of Change of 

the Control Volume “ the Control Volume + Energy within the 



Control Volume 



mh + 



( 





or, upon simplifying the above equation. 
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In equation (A.l), the first term is the change in enthalpy flux 
across the control volume, the second term represents the gradient of 
the energy transport by gas conduction, the third term is the convec- 
tive heat transfer between the gas and the matrix, and the fourth term 
is the rate of change of energy within the control volume. 



matrix Energy Equation; 



Energy flux into _ Energy flux out of Rate of change of 

the Control Volume “ the Control Volume Energy within the 

Control Volume. 



*^Sv, A + hrpAm 

*“ £ x 1 1 



T - T dx s 



^ A H + k » A £|£| & 



a 



+ Me — dx 
» it 



or, upon simplifying. 



ICja A + bjAj \ - T 



- Kc 



m 



dT 



(A. 2) 



In equation (A *2) the first term is the gradient of the energy 
transport by conduction along the matrix, the second term is the con- 
vective heat transfer from the gas to the matrix, and the term on the 
right is the rate of change of energy of the matrix. 

Equations (A.l) and (A .2) are energy balances for the gas and 
matrix, respectively. Subscripts g and m refer to the gas and the 
matrix, respectively; h represents enthalpy and u represents internal 
energy; h^ is the heat transfer coefficient (which has not been assumed 
constant) and k is thermal conductivity; A is the total cross-sectional 
area of the regenerator and A 0 is the cross-sectional area open to 
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flow (a o z € A); is the heat transfer surface area of the regenerator 
per unit length (ft^/ft)* and M is the matrix mass per unit length (ib^/ft)* 
The gas continuity equation will now be developed. 

Gas Continuity : 



Mass flow into the 
Control Volume 



Mass flow out of the 
Control Volume 



Rate of storage of 
mass within the 
Control Volume 



0 


. * \ 


+ JL 


m - 


m + uJE dx I 




dx 1 


b t 



<? A dx 

go 



} 




dra d ?cr 

dx A 0 bt 



(A .3) 



Equation (A. 3 ) is the continuity equation and it is perfectly general 
(within the assumptions listed previously), as it applies to either 
steady or unsteady flow. With the use of the perfect gas law. 



9 Z - , (AJi) 

S P.Ig 

equation (A . 3 ) becomes 

b m _ _ A 0 ^ 
b x R 3 1 

Equations (A 0 l), (A.2), and (A 0 3) along with appropriate boundary con- 
ditions are the governing differential equations for any type of regen- 
erator operation. Equation (A.l) will now be simplified by using 
equations (A 0 Ii ) and (A„5). 

For a perfect gas, one may write 



‘•2 



(A.S) 
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(Ac 6) 



h - h r 



* c P T g 



and 



u = h-pv=h-RT = h * c T , 
^ g o v g * 



where h Q is the gas enthalpy at some reference temperature • 

The first term in equation (A.l) becomes, with the use of equa- 
tions (A. 5) and (A »6): 



dx 







_ AqIIq 6 

~R dt 






a 



(A.7) 



The fourth term in equation (A.l) becomes, with the use of equa- 
tions (A .U) and (A .6): 



* i- f <? u\ zt£-° -I?-) ^ . (A.8) 

0 dt \ S ) R dt \ T g J R dt 



The substitution of equations (A.7) and (A.8) into equation (A.l) 
yields : 




- k A 

E 7i? 



hrnA r-» 



« T 



c- yA-o 

R 



6 P 
<H 



- 0 



(A.9) 



Equation (A.9) is the final form of the energy equation for the gas. 

Equations (A.9), (A.2), and (A.5), along with the appropriate 
boundary conditions, completely describe the performance of thermal 
regenerators consistent with the assumptions listed on page 57. 

The assumption will now be made that terms involving thermal 
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conductivity are negligible in comparison to the other terms in (A 0 9) 
and (A 0 2)* That this assumption is consistent with some practical 
cases is shown by taking some typical valuesr 

k g = 0.07 ETU/hr-ft-°F, A Q z 0 o 001$k ft 2 , (T„ - T) = 2°F, 

A t = U.1j 3 ft 2 /ft, hm z 235 BTU/hr-ft 2 -°F, 6 2? g _ T W“ T o - i5o°F/ft 2 . 

b x 2 L 2 



Then: 

A 2 T 

Vo 6 g " 0.0201* BTU/hr-ft, and 

^ " 

hpAy (T g - T) = 2080 BTU/hr-ft o 

Thus, the gas conduction term is approximately 0*001$ of *che con- 
vective heat transfer term of equation (A *9). 

In the matrix energy equation, (A«2), with k m - 20 BTU/hr-ft-°F 
and A = 0.005 ft 2 : 

A 2 T 

k A — 15 BTU/hr-ft* 
m ^2 

Therefore the conduction term for the matrix is less than 1$ of the con- 
vective heat transfer term. Then neglecting axial conduction, the final 
basic differential equations become: 



c 



P 



_6_ 

x 




h/pAjp 



g 



c v^o b P 

IT ft 



o 



(A *10) 







= 

m 



dj? 

At 



(A*ll) 
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p 




(A .5) 



b x r n t 

' I 



g 



In equation (A.10), the first tern is the gradient of the enthalpy 
flux along the regenerator; the second terra is the energy flux due to 
convective heat transfer ; the third terra is the work of compression. 
Equation (A .11) means that the increase in internal energy of the 
matrix equals the energy given up by the gas by convective heat trans- 
fer. Equation (A .£) is the gas continuity equation. 

Boundary conditions for the above differential equations are: 



Simplif ic at ion of the Differential Equations 

Equations (A. 10) 3 (A.H), and (A 0 £) are three equations in three 



unknowns: T, T^, and m. Thus, if p as a function of t is known, then 

solutions consistent with the boundary conditions, equations (A .12), 
will exist. Unfortunately, equations (A .10), (A .11), and (A .5) are a 
system of differential equations for which no general solution is known. 
The method of approach, then, will be to determine solutions which 
satisfy certain conditions over a complete cycle, but which are not 
exact in every detail at each instant within the cycle 0 

A cyclic integral of equation (A .11), while holding x constant, 
yields : 




(A. 12) 



at x Z L: T 




© 
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r 

© hAm 
; * 



(A.13) 



‘g 



QC 




ST 

St 



dt = 0, 



The above follows from the fact that once steady, cyclic conditions 
have been established, the matrix temperature at a given point in the 
regenerator oscillates around a value of the temperature which is a 
function of x only* From this point on, h will be used for the heat 
transfer coefficient rather than h^« 

A cyclic integral of equation (A *10) with x constant yields: 




dt v fzjo dt = 0* 
R ) b t 



(A.1U) 



But, 



b t 



dt - 0 by definition, since p is a function of t only 



(no pressure drop in the regenerator)* Then by using equation (A.13), 
equation (A.lU) becomes: 



h h T g dt °- 



u.i;) 



« 

Since nT^ will be continuous in x, the order of differentiation and 

o 

integration can be interchanged and equation (A.l£) becomes: 



6 x 







or, 

mTgj dt - C 0 (for all x). (A«l6) 

For the 2-Part Model, the pressurization part of the cycle is 
denoted by the subscript c and the expansion part of the cycle by the 
subscript e* The mass flow rate at a given point within the regenerator 
will be cyclic, passing through zero at reversal. The mass flow versus 




-6 



o- 



time might appear as: 




Since there can be no net mass storage within the regenerator once 
steady, cyclic conditions prevail. 



f 



mdt 



0 



(A .17) 



and the positive area must exactly cancel the negative area in the 
above diagram* Similarly, will be cyclic with time, but in some 

o 

unspecified manner and (f T g dt will not, in general, equal zero. 

o 

Then, the product mT g might have a shape as follows: 





(nT g ) e 

where the mT g dt - C Q from equation (A.16). This means that the 
positive and negative areas in the above sketch do not, in general, 

o 

cancel each other. The assumption will now be made that mT g versus 
time can be represented by the dashed lines in the above figure o 
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Then? 



mT g dt = ) c rn c T gc dt c * J>Vge dt e = c 0 



- m c T gc r c + m e T ge r e = c 0 



(A.18) 



where / T is the blow tine during compression and / T Q is the blow tine 
during expansion,, 

Furthermore, if I' - T" (compression period equals expansion 

c © 

period ), equation (A .18) can be written: 



o o 

m T + m T 
c gc e ge 



(A .19) 



where 



°1 = c o/T c = C 0 /^ £ 



From equation (A 0 17) s 



mdt Z 



m c dt + \ m e dt " 0 



md 5 



|c m C dt = ) e ra e dt 



(A*20) 



With the assumption that the mass flow rate can be represented by an 
average flow rate which is constant with time over each half-cycle, 
equation (A„20) becomes: 



-c T c = - 



m e T e 
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(A. 21) 



and again letting 'T z T 

c e 



m c = 



Equation (A*20) and, hence, equation (A. a 2l) holds at any given x 
along the regenerator. From here on, from the result of equation (A .21), 

o 

m id.ll be used without subscript to denote flow in either the compression 

a 

or expansion half-cycle and n is a function of x only* 

With the assumption that 



mT = (m)(T) , 

equation (A e 19) along with (A 0 21) becomes 



™ ( T gc " T ge) = C 1 



(A *22) 



Equation (A 0 22) is a very important and useful result. There are no 
restrictions on C-, * It can be positive, negative, or zero, depending 
upon whether T is greater than, less than, or equal to To 
The usual heat transfer correlation is assumed, 

h = V 1 > 



where n is a positive exponent whose value depends on the particular 
regenerator matrix. In general 0*£ — n — 0.7. Hence, h is either 
constant or a function of m only. Using equation (A .21), h xdll be a 
function of x only and will be written as h. Then from equation 
(A *13): 



h A (T - 
T g 



1 



T) dt - h A r 



(T g T) 



- 0 



since h is independent of t, or. 
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(Ao23) 




A plot of (T g - T) versus t might appear as: 




From equation (A .23), the positive area on the above sketch must equal 
the negative area. As before, assuming (T - f) can be represented by 

o 

some average value and letting T c r 'T' , equation (A.23) becomes: 

(T g - T) c = (T g - T) e (A.2U) 

Since the blow-time, 'X , is usually short, a reasonable assump- 
tion to make is that the temperature of the matrix at a given x will 
vary very little during a half-cycle and, in fact, T can be assumed 
constant with time. In that case, T - T (x). The validity of this 
assumption can be demonstrated by the following example. From equa- 
tion (A.ll): 

h A T (Tj, - T) = Kc m £2 ( Aoll ) 

For a regenerator 12 ,; long and 1" diameter packed with lead shot of 

0.050" diameter, Mc^L = .068 BTU/°F. Taking h r 235 BTU/hr-ft 2 -°F, 

A^L z i;.U3 ft 2 , and assuming (T - T) 2°F: 

bA„ (T - T) i 2080 BTO /hr I 0.068 *2 

1 ° d t 3 
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fivh a blow time of 0.1 second,AT — 0.85 °F 0 During a complete cycle, 
at a given x, the matrix temperature would oscillate around a m^an tem- 
perature with an amplitude on the order of 0.85 °F. Thus, the assumption 
is made that T is independent of time and is a function of x only. 

By applying the above assumption to equation (A 0 2li)> one is lead 
to the important result that: 



At first glance equation (A .25) seems obvious, but a careful examina- 
tion of the preceeding discussion will show that two key assumptions 
were necessary in deriving equation (A. 25). The assumptions wore 



The ratio of the first term in the brackets of equation (A. 26) to the 
second term is: 




(A.25) 



(l) the mass flow rate, m, is a function of x only, and (2) T is a 
function of x only. 



By carrying out the time-differentiation in the continuity 
equation (A .5), one obtains: 




o 



1 dp P dT AP T T 

rp 8p /-w p-y 

T dt ? 2 dt p AT„ AT 

g g S g 




With 



AT -AT 
g 



0.85 °F , the above ratio becomes: 
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U30 to 625 






AT 



g 



1.2T z 

g 



Thus, the second term in the brackets of equation (A 0 26) is less than 
1% of the first term and, consequently, the second term will be 
neglected. Equation (A*26) can now be written as: 



dm A Q dp 
dx " RT dt 

o 



(a *27) 



Total derivatives are used in equations (A 0 27) since p is a function cf 

a 

t only and m is considered to be a function of x only. Furthermore, T 

o 

trf.ll be replaced by T in equation (A.27) since (T - T) is small 

(2° F to 5° F). This assumption also enables continuity to be satisfied 

0 

over a cycle. Since m has already been assumed to be equal and opposite 

for the compression and expansion half-cycles, it is reasonable to expect 

Sri to be the same for both half -cycles. Then the continuity equation 
ccc 

is written as: 



dm _ A 0 dp 
dx RT dt 



(A 0 28) 



It is now assumed that S r constant. Thus, pressure is 

QU 

assumed to be a sawtooth function of time as shown in the sketch below, 
where the zero-axis is any arbitrary pressure level: 
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+ constant also maintains consistency with the previous 



Letting ^5 r 

* 

assumptions that m and T are functions of x only* 

It is now possible to present the final formulation of the regen- 
erator problem. From Figure XII, the positive x direction for the 
regenerator is from bottom to top. Mass flow in this direction is posi- 
tive. During pressurization the flow is assumed to be from top to 
bottom, or, in the negative x direction. The gas energy equation which 
was derived for positive m with positive , is written for the two 
half-cycles: 



-c 






-T A gc 



< 

i 

o 


dp 


1 R 


dt 



= 0 . 



(A .2 9) 



c pl_l^ T ge| + M T 


> - t I 


- CyA 0 


dp 


dx ' / 




R 


dt 



I 0 . 



(A. 30) 



In equations (A„29) and (A o 30) 3 



dp 

dt 



is written as an absolute value 



with the plus or minus sign placed in front, depending upon whether 
pressurization or expansion is occurring in the regenerator. The con- 
tinuity relation, equation (A. 28), is written: 



dm 

dx 



RT 



£P 

dt 



(A .31) 



For convenience, average value bars have been omitted. Equation (A. 31) 

o 

holds for both half-cycles, since in the compression half-cycle m is 
negative, but £P is positive, whereas in the expansion half-cycle 

Q X 

o 

m is oositive, but is negative. 

dt 

The final formulation of the problem consists, then, of equations 
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(A 0 29), (A.30), snd (A. 3l) along with the results previously attained: 



m 



T 

£ c 




(A .22) 




2T 



(A.25) 



Boundary conditions for the above set of equations are: 

at x : 0: T z T 0 

o » 

m - m Q 

at x = L: T = % . (A 0 32) 

Equations (A 0 29), (A.30), (A*3l), (A„22), and (A. 23) are five 

* 

equations in five unknowns: m, T, Tg C , T 1 , and However, there 

are only three boundary conditions. Therefore, it is necessary to 
eliminate the gas temperatures from the above equations. 

If the differentiation indicated in equations (A. 29) and (A.30) 
is carried out there results: 



~c T 



l £C 



dm 



m — £X + Tgc — + hAj T gc - T 

\ dx ox 1 1 



°-'Ac 

R 



dt 



- 0 



(A .2 9a) 






m + T ge — + M T I T ge “ T j “ ^ 

dx dx / \ ' R 



dt 



~ 0 



(A. 30a) 



Now if equation (A. 30a) is subtracted from equation (A*29a) the result 
is: 



-C v 



m + £Ig® 

dx dx 
2®v^o 



1 dm 

* T s e 



dp 

dt 



+ hJl T T gc " T ge 



= 0 



(A.33) 
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Now by substituting equations (A, 22), (A„23), end (A.31), plus the 
relationship R : Cp - equation (A. 33) reduces to: 

m ^ _ b^'T 0 ! + ^0 

dx 2c n c 
P P 



dp 

dt 



1 0 



(A .31) 



which involves only the matrix temperature, T. 

The following new variables are now introduced: 



© : 



iA 0 



m 

w 




V dp 

T m R dt 
0 o 



X 

L 



(a.35) 



Also the heat transfer coefficient is now assumed to be correlated by 

h = Km 1 , (A. 36 ) 

where K is a constant with the units BTU/hr-f t -°F. 

With the aid of the new variables, equation (A. 36 )* and the 
definition z c p , equation (A. 3 I 1 ) can be written: 



m dw y 



1-n 



m 



C 



where 



C = 



2m_.V 


dp 




0 


dt 





The continuity relation, equation (A 0 3l)> becomes 

dm - 1 

dw “ © 



(A. 37) 



(A.38) 



(A. 39) 
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Boundary conditions for equations (A. 37) and (A«3 9) are: 



at w : 0: si: 1 

0=1 

: © = © w (A.IjO) 

Equations (A *37) and (A. 39), along with the boundary conditions 
(A o l|0), are sufficient to determine m, T, and the constant C. Actually, 
since everything in the definition of C is known except C^, determining 
C is essentially determining C^. 

Thus the problem for the regenerator is reduced to one of solving 
a set of two ordinary differential equations subject to the boundary 
conditions mentioned. The method of solution will be discussed later 
in section A. 3. 



at w z v 



max 



V 

T 0 m 0 K 



do 

dt 
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THE GOVERNING DIFFERENTIAL EQUATIONS FOR 



REA»S 2 -PART MODEL AS APPLIED TO A " pulse tube " 

A pulse tube, as long as its heat transfer area per unit total 
volume is fairly large, can be considered to be an inefficient regen- 
erator. Thus, a similar development to that leading to the final 
differential equations for the regenerator could be presented for the 
pulse tube. However, such a development would be repetitious. The 
procedure that will be followed will be to accept the assumptions 
listed on page 57 and start with equations (A.9), (A.5), and (A .2), 
which are the gas energy equation, the gas continuity equation, and 
the matrix energy equation respectively. Appropriate assumptions 
regarding the relative magnitude of the axial conductivity will then 
be applied. Next the cyclic conditions will be presented. At this 
point the relative magnitude of the terms of the continuity equation 
will be examined. The final formulation of the pulse tube problem 
will then be possible. The pulse tube differential equations will 
evolve in much the same form as equations (A. 37) and (A, 39). It 
will be of some assistance to present here a schematic for the 
2-Part Model for a pulse tube. This is shown in Figure XIII. Note 
that x is measured positive from the top down, which is directly 
opposite from the case of the regenerator. 

Equations (A .9), (A.5), and (A .2) in the form in which they 
were previously developed are as follows: 



(A.9) 



c p i- 



mT 



- Vo £? ♦ “i ( T * 



. t .^*2:0 

I R dt 



d m 



ax R dt 



(Ao5) 



FIGURE Xin 



SCHEMATIC OF ITffi SYSTEM TO BE ANAITZFD FOR THE " pulse tube " 



A * 

i 



Control f 
Volume 



max 



T dx 



P . 
min 



A 



+ Vt 






T - T = Me 



m 



ajr 

at 



(A .2) 



As in the case of the regenerator, the assumption will now be made 
that the terms involving thermal conductivity are negligible in compari- 
son with the other terms in equations (A. 9) and (A.2). Some typical 
values which substantiate this assumption are as follows: 
k r 0.07 BTU/hr-ft-°F, A = 0.0019 ft 2 , (T - T) = 2°F, 

6 ^ o 
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A t = 0.155 h T = 6.18 ETU/hr-ft-F, 

. 2rn rn o 

o ^ A H “ x o _ 2 

— - 120°F/ft . 

8 x 2 L 2 



Then: 

2 

lc g A o ^ Tfr : 0.016 BTU/hr-ft, and 

i x 2 

V^T ( T g ~ T ) = 1.92 BTU/hr-ft. 



Thus, the gas conduction term is less than 1 % of the convective heat 
transfer term of equation (A *9)« In the matrix energy equation (A. 2), 
with 10 BTU/hr-ft-°F and A = 0.000225 ft 2 : 

\ 2ip 

~ 0.270 BTU/hr-ft 

byr 

the conduction term is approximately ll$ of the convective term. If the 
device were operated with a (T^ - j^) of 190°F, this percentage would 
increase to approximately 22/. Obviously this assumption is not as good 
as in the regenerator case, but without it the problem would become 
exceedingly complex. Therefore equations (A.9) and (A .2) become 
respectively: 



_ A I mT„| ♦ hAm T - t| * CyA ° : 0 
E ) T \ g I R d 



and 

h T^T [ T g 



\ 

T 



r Me 

m 



dr 

6 1 



(A .10) 



(A. 11) 



The cyclic conditions for the pulse tube are identical with those 
for the regenerator. Thus, equation (A. 22) applies equally as well to 
the pulse tube as to the regenerator. Again the heat transfer 
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correlation is assumed to be: 




For the pulse tube, however, laminar flow is assumed, /t the moment it 
is not intuitively obvious that this is a valid assumption. However, 
experimental data indicates that laminar flow is probably the closest 
approximation that can be made for the pulse tube. Thus, n = 0, and 
the heat transfer correlation for the pulse tube becomes: 



Now, proceeding with an identical argument as that used in the case 
of the regenerator, one obtains equation (A.210: 



Here the assumption is made that the temperature of the pulse 
tube wall at a given value of x will vary little during a half -cycle. 
Therefore T can be assumed to be constant with time, and T = T (x) only. 
The validity of this assumption can be shown by the following example. 
From equation (A .11): 



For a stainless steel pulse tube 12” long with O.D. - 0.625" and I.D* = 
0.591", Mc^L = 0.011 BTU/°F, and taking h = 6.18 BTU/hr-ft 2 -°F, 

AjL = 0.155 ft 2 , and assuming (T g - T) 2°F; 



h = constant 




(A.21) 



hAm( T - T) r Me 41 
1 8 m A*. 



(A .11) 



hA_,(T - T) = 1.92 BTU/hr = 0.011 4J , 
I g o t 



and 
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Thus, with a blow-time of 0.1 second, A T — 0.005°F. During a com- 
plete cycle, at a given x, the pulse tube wall temperature would 
oscillate around a mean temperature with an amplitude of approximately 
0.005°F. 

Thus, the assumption that T is independent of time and is only a 
function of x, is a valid one. 

By applying this assumption to equation (A.2lj), one obtains, once 
again, the result: 



T gc +T ge=' 2T • 



(A.2S) 



By carrying out the time differentiation of the continuity equa- 
tion (A.5), one obtains: 



• 

d m 

mm 




'1 


a p 


p 


d X 


■r 


T 


at 


t 7 



b t 



(A. 26) 



The ratio of the first term in the brackets of equation (A .26) to the 
second term is: 

1 <3p / V A? T 

a ZL — — .ia 

t„ at/ T g at p AT g 



From equation (A,22)r 

AT - ^1 
m 

From equation (A. 38 ): 



(a. la) 



C, r 2C -i~ mo 



V 



-1 KAipL 



dp 

dt 



(A.U2) 



By substituting equation (A.lj2) into equation (A.Ul)> ope obtains: 
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V 



dp 

dt 



(A.U3) 



m. 



A Tg = 2C _£_ 

$ -1 m KA t L 



Thus, 

AP Tg 1 #-l m KArpL At 

p AT g 2C y m Q V p Tg 

2 o 

Now using typical values of s 1*67, KArjL/V - 1*1.6 BTU/ in-hr-ft - F, 

C - ,2, m/m 0 = 3> At- .27 sec, p ; 75 psia, T g = 3$>0°F, the above 
ratio becomes: 



A J- Is - 3 

P AT g 

For the regenerator this same ratio was U30 to 62 $ , and the second term 
was neglected. Rea’s theory also neglects this term for the pulse tube. 
However, the above example indicates that the second term cannot be 
neglected in the case of the pulse tube. Equation (A ,26), then, remains 
as: 

dm _ A 0 / 1 dp p d^ 

dx R \T dt T 2 dt 

\ s s 



(AM) 



Total derivatives are used because p is assumed to be a function of 

t only, T (at a given x) is considered to be a function of t only, and 
& 

m is considered to be a function of x only. Also, since (T - T) is 

b 

small ( 2°?. to 5°F.), T will be replaced by T in equation (A.l*3)» Then 

© 

the continuity equation can be written as: 

dm A q I 1 dp p 

dx R T dt T 2 



dT_ \ , 

— s . (AM) 

dt / 
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Now the assumption is made that dl g can be approximated by 

dt 

o Once again is assumed to be constant for both half-cycles, 

£ tr dt 

Now by substituting equation (A.U3) into equation (A.lt£) and using the 

fact that + constant, equation (A.lj£) becomes: 

at - 



dm _ A q 



dx R 
or, simplifying, 



dp 

dt 



- 2C 



_L] L 1 W 

2T-1/ AtWlm 



— 

KA^/ 



dp 

dt 



(A.U6) 



• 

dm _ A 0 


£p 


H 


dx "* RT 


dt 


where 






B = 2C j 


f j_' 


1- - 


1 


itf-u 


! At T 



m- 



in 



ka t l 



(A.lt7) 



(A.ii8) 



In equations (A.U6) and (A t Li7 ) average value bars have been omitted for 
convenience. 

As in the case of the regenerator, the gas energy equation can be 
written for the two half -cycles: 



—C r 



ml 



ax 



gc 



+ hAr 



l gc 



T 


+ 


dp 




R 


dt 



0 



(A.29) 







+ hAfp 





CyA 0 dp 
R dt 



I 0 



(A.30) 



The final formulation of the pulse tube problem consists, then, 
of equations (A.U7), (A. 29), (A. 30) along with the results: 
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m 



9 



(A. 22 ) 





gc 



*T ge r 2T 



(A.2S) 



Boundary conditions for the above set of equations are, as for the 
regenerator: 

at x : 0: T = T Q (A.U9) 

at x : L: T = T c 

Note that the condition at x : L is different for the pulse tube than 
for the regenerator only in the fact that the cold end of the pulse 
tube is at x Z L whereas the warm end of the regenerator was at 
x = L. Thus the subscript C stands for "cold” in the pulse tube bound- 
ary condition, whereas the subscript W stood for "warm" in the regen- 
erator boundary condition,, 

As in the case for the regenerator, equations (A. 29), (A„30), 
(A*U7), (A«22), and (A, 25) are five equations in five unknowns: 
m, T, Tg C , Tg fi , and , Since there are only three boundary conditions, 
it is necessary to eliminate the gas temperatures from the above equa- 
tions • 

If equation (A. 30) is subtracted from equation (A. 29) and if 
equations (A.22), (A, 25), and (A'.U7) are used to eliminate the gas 
temperatures, the resulting equation involves only the temperature of 
the pulse tube wall: 



ro dT _ hA^C x + A 0 


dp 


i 


[l-B ' 


1* 


dx 2Cpm Cp 


dt 




y 


-l 



(A.50) 
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Nov if B is replaced by equation (A.I 48 ), equation (A o 50) becomes: 



• dT 

m 



hA^C^. A o 



dx 2Cpi5i 



dt 



1-2 C M P 1 K 



V 



\\ 



U-i/ At t 



m 



,KA t L/ 



J-i 



* -1 



- 0 • 

(A.5'1) 



The variables ©, m, and w, equation (A ,35), plus the heat trans- 
fer correlation, h = K, are now introduced into equation (A 0 5l)> 
obtaining: 



a 2 :c 

dw \ 2r 1 



1 

© 



where 



<J> = 2 



y-i, 



v P_ i_ 

^ KA.ipXi| Z\ t 



(A .52) 



(A. 53) 



A t is £t „ for the entire cycle, and p is p for the cycle which is 
avg 7 avg 

taken as: 



p - p . + ^max “ ^min 

avg " min 

2 

The continuity equation (A. 1*7) becomes in this case: 



(A.5U) 



dm _ 1 
dw “ Q 



1 t£ 

x ‘©m 



(A. 55) 



The boundary conditions for equations (A<,52) and (A 0 5U) are: 



at w : 0: 



m Z 1 

Qz 1 



(A .56) 



at w - w, 



_ V 



max 



T m R 
0 o 



8 p 

dt 



■ e-.e_ 
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Thus, equations (A, 52) and (A, 55) together with the boundary 
conditions (A.56) describe the problem for the pulse tube. 






A.3 



METHOD OF SOLUTION FOR THE DIFFERENTIAL EQUATIONS 



Regenerator 

If equation (A. 37) is solved for , 

dw 



d© 

dw 




2-n 



m 



Also, 



dm - 1 

dw © 



(A .57) 



(A. 39) 



The simplest method of solution for the set of equations (A. 57) 

and (A *39) is to solve them numerically with computer aid. The 

simplest approach is to treat C as a parameter and assign a range of 

values to it. The problem then becomes an initial value problem, since 

everything on the right side of equations (A. 57) and (A. 39) at w Z 0 is 

known. The independent variable w can then be run out to sufficiently 

large values to include any desired w^^. 0^ then becomes an unknown 

quantity, but can be determined from a plot of 0 versus w for a given 

value of C. On the other hand, for a given value of © TT and w , the 

3 b W max* 

corresponding value of C can be read from the same plot. 

Figure XIV shows a hypothetical set of solutions corresponding . 
to a family of C’s. From the figure, for C = Cg, © = at 
w = w . The entire regenerator -temperature distribution is deter- 
mined by following the curve for C S CL from w = 0 to w = v8 In 

£ lucLX 

like manner one can generate curves of m versus w from equation 
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(A. 39). In this case the mass flow distribution can be determined by 

following the appropriate C curve from w : o to v : w . In addition 

max 

to specifying a range of C’s, the regenerator matrix must be known as 
well as the working gas, so that values of n and % can be assigned in 
equation (A.57)* 



FIGURE XIV 

ILLUSTRATION OF COMPUTER SOLUTIONS FOR REGENERATOR 




- 8 $- 



" pulse tube " 



If equation (A. £2) is solved for 

dw 



dQ _ 
dw 



C 



i 

1 + 




nr 




Also 



dm _ 1 $C 

dw ** © \ ©m / 



(A .58) 



(A.55) 



Again, the simplest method of solution of these equations is by 
computer. Similar to the case for the regenerator, the simplest approach 
is to treat ^ and C as parameters. Thus, for a given (j) z 
C is assigned a range of values j then a new value is chosen for and 
C is assigned another range of values. The result is a family of curves 
of C for a given value of <£ , as is shown in Figure XV. The curves are 
used in exactly the same manner as that described for the regenerator* 

To use the curves one must be able to determine in advance. From 
equation (A.5>3) it is apparent that to determine ^ , the following 
information must be known: 

(1) the working gas, 

(2) the geometry of the pulse tube, 

(3) the cycling rate (See Chapter IV), 

(ii) the temperature of the cooling medium, 

(5) the average pressure at which the device will be operated 
(defined by equation (A.5ii)). 



Details of the computer program for the regenerator equations can 
be found in (ii). Details of the computer program for the pulse tube 
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equations can be found in Appendix B. Computer solutions for the regen- 
erator and pulse tube can be found in Appendix C* 



FIGURE XV 

ILLUSTRATION CF COMPUTER SOLUTIONS FOR " pulse tube " 





w 



— > 

w 



APPENDIX B 



COMPUTER PROGRAM FOR THE SOLUTION OF THE "pulse tube " EQUATIONS 
cncral 

Rea (1|) wrote a program in Fortran to solve the regenerator equa- 
tions (2.2 9) and (2.30). A similar program was written to solve the 
pulso tube equations (2.33) and (2.3h). Runge's recurrence formula (6) 
was used in the stepwise integration. 

The program was written for X z 1.67. FUNF(B,R) is equation (2.33 )> 
where B corresponds to the mass flow, m, and R corresponds to the dimen- 
sionless temperature ratio, © . The other symbols used in the program 
are: 

PAR = C 
SPAR = $ 

TEMP = © = T/T q 

Y - m = m/m 0 
SLOPE = d©/dw 

D - m* in Runge recurrence formula 

E = ©* in Runge recurrence formula 

H = integration step size = 0.01 

K = number of integration steps 

XK = number of values of C for each value of $ 

KKK = number of values of <£> 

The first data card contains the initial values of © and m, the 
step size, the number of integration steps, the number of values of C 
for each value of , and the number of values of (£) . The remaining 
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data cards contain the values of <§> and C for which results are desired. 



Computer Program 



* XER 

* LABEL 

* LISTS 

* FORTRAN 

FITMF ( B , R )= ( PAR* (l . +SPAR/R ) - ( 0 ,I>*B ) )/( B*B ) 

DIMENSION A ( ),C( ),TEMP( ),Y( ),SLOPE( ),W( ) 

READ 10,TEMP(l),Y(l),H,K,KK,KKK 
DO 2^ L=1,KKK 
READ 12, A (L) 

SPARzA(L) 

3 DO 25 1=1, KK 

READ 12,C(I) 

PARsC(I) 

6 DO 20 J=1,K 

D=Y(j)+(H*(l./TEMP(j))*(l.-((SPAR*PAR)/(TEMP(J)*Y(j))))/2.) 
E=TEMP(J ) +( H-*FUNF ( Y( J ) , TEMP ( J ) )/2 . ) 

TEMP ( J+l ) =TEMP ( J )+ ( HftFUNF ( D, E ) ) 

Y( J+l )= Y( J ) +(H* (1 o/E)*(l . - ( ( SPAR*PAR )/( E*D) ) ) ) 

¥ ( J+l )=H*FLOATF ( J ) 

20 SLOPE ( J +1 )=FUNF (Y( J+l ) , TEMP ( J+l ) ) 

PRINT 13, H, PAR, SPAR 

25 PRINT 15, (W(J+1 ),TEMP(J+l),Y(J+l),SL0PE(J+l),J=20,K,20) 

CALL EXIT 

10 FORMAT (3F12.6„IU,I2,I2) 

12 FORMAT (1F12.6) 

13 FORMAT (3F12.6) ' 

15 FORMAT (iiFl5o5) 

END 

* DATA 

1.0 1.0 .01 600 710 
(OTHER DATA FOLLOWS) 
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APPENDIX C 



THEORETICAL AXIAL TEMPERATURE 
AND AXIAL MASS PLOT DISTRIBUTIONS , 
2-PART MODEL. 
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APPENDIX D 



DESCRIPTIO N OF THE TEST APPARATUS 

The test apparatus consisted of an inlet/exhaust manifold. Pulse 
Tube, heater, heat exchanger, and vacuum chamber. A schematic of the 
apparatus is shown in Figure IV of Chapter IV. 

The inlet/exhaust manifold is shown in Figure XXIX. The manifold 
was constructed of standard copper-tubing, using both soft soldered 
and flared fitting connections,. Four Asco solenoid-valves were used to 
control the flow, two on the supply side and two on the exhaust side of 
the manifold. These valves had 3/32" orifice. The solenoid valves were 
cycled by a cam operated micro-switch arrangement. The cams were belt- 
driven by a small, electric motor through a four-step pulley which 
allowed cycling rates of 22 ,U, Uu7, 82.3, and 158.0 reversals per minute. 
Two 3/8" Koke needle-valves, one upstream of the inlet solenoid-valves 
and one downstream of the exit solenoid- valves, provided fine control 
of the gas flow rate. The check valve/manual valve loops allowed for 
reading maximum and minimum pressures without any pressure gauge 
fluctuation. The maximum pressure was read on a !;"-Marsh, 0-300 psi 
gauge. The minimum pressure was read on a 3''“Marsh, 0-100 psi gauge. 

The copper- tubing of the manifold mated to a 1" O.D. by 2^-"-long, hollow, 
brass cylinder. The lower end of this cylinder was machined to allow 
an "0" ring to be inserted. This end of the cylinder was fitted into 
the open-end of the regenerator. Several short lengths of tubing were 
fitted inside the brass cylinder to act as flow smoothers. 

The Pulse Tube drawn to scale is shown in Figure XXX. It consisted 
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FIGURE XXX 



SCALE DRAWING OF PULSE TUBE 
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SCALES = 1,0" 



of a regenerator, a void tube (pulse tube), and a short length of connect- 
ing tubing. The regenerator was constructed from a 12” -Ion g, stainless- 
steel tube, 1" O.D. by 0.020" -wall- thickness packed with 0 o 650" -diameter 
lead shot. The pulse tube consisted of a 12"-long, stainless-steel tube, 
5/8" O.D. by 0.0l6"-wall-thiclcness. The regenerator and pul.se tube were 
connected at the bottom by a -v" O.Do by 0 0 02 0 " -wall - 1 hi c kr.e s s , stainless- 
steel tube. On the pulse tube end of the connecting tubing was a 1" I.D, 
by l"-long, stainless-steel cylinder, which was filled with fine-mesh, 
bronse screens to serve as a flow smoothing device. Several layers of 
this screening were also placed at the bottom of the regenerator for the 
same purpose 0 The stainless-steel components were heli-arc welded as 
shown in Figure XXX. The top ends of the regenerator and pulse tube 
were soft soldered 3" apart to a brass, vacuum- jacket, cover-plate, which 
was 9" 0 o D. by 3/8 !, -thick 0 Eleven copper-constantan thermocouples were 
installed along the regenerator and pulse tube at locations indicated in 
Figure XXX, All the thermocouples were led out of the vacuum chamber 
through two l"-diameter Latronic ceramic seals, which were soldered into 
pre-drilled holes in the cover plate. The thermocouple outputs were 
recorded by a Leeds and llorthrup multipoint recorder. The pressure- 
time history within the Pulse Tube was measured by diaphragm-type 
pressure transducers mounted at the top of the regenerator and pulse 
tube. An oscilloscope was used to indicate the output of the pressure 
transducers. 

The heater was constructed as follows: the jj" stainless-steel 

tubing connecting tho regenerator and pulse tube was coated with a thin 
layer of "Sauereisen Insa-Lute Adhesive Cement", a liquid-porcelain type 
of material. After the cement had hardened, approximately two maters of 
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a high resistance wire were wrapped on top of it* The heater wire was 
connected to two 32-gauge, copper, lead-in wires. Then the entire 
heater assembly was coated again with a second layer of cement. The 
lead-in wires were lead through a Latronic ceramic seal installed in 
the vacuum-jacket cover plate. The total resistance of the heater cir- 
cuit was 121 ohms. The heater circuit external to the vacuum jacket 
consisted of a O-lljO volt variac for fine voltage control 0 A milli- 
ampmeter, voltmeter, and wattmeter were installed in the external cir- 
cuit for power measurements. 

The heat exchanger at the top end of the pulse tube was constructed 
as shown in Figure XXXI. Thin-wall, copper-nickel tubing l/8" 0,0. was 
used to promote good heat transfer. The working gas was funneled into 
the tubes, which were approximately in length and which were dead- 
ended in a brass tube sheet approximately 1'* in diameter by -| t, -thick* 

The entire tube assembly was surrounded by a water jacket, which was 
attached to the cover plate. Cooling water was pumped through the water 
jacket at the rate of 0,5 gallons per minute. 
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FIGURE XXXI 



HEAT EXCHANGER 
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